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^ ■ Abstract 

h : 

Almost all research on superintegrable potentials concerns spaces of constant 
curvature. In this paper we find by exhaustive calculation, all superintegrable po- 
tentials in the four Darboux spaces of revolution that have at least two integrals 
of motion quadratic in the momenta, in addition to the Hamiltonian. These are 
two-dimensional spaces of nonconstant curvature. It turns out that all of these po- 
tentials are equivalent to superintegrable potentials in complex Euclidean 2-space or 
on the complex 2-sphere, via "coupling constant metamorphosis" (or equivalently, 
via Stackel multiplier transformations). We present tables of the results. 
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1 Introduction 

In a previous article £Q we have studied superintegrability in a two-dimensional space 
of nonconstant curvature, in particular one of the so called Darboux spaces, given by 
Koenigs [2] . In this article we study the remaining three spaces of nonconstant curvature 
from the point of view of superintegrability. This involves the addition of a potential to 
each of the spaces given by Koenigs. We recall that classical superintegrability relating to 
a Hamiltonian H(x\, . . . , x n ,pi, . . . ,p n ) = H(x,p) implies the existence of 2n — 1 globally 
defined constants of the motion. For the purposes of this article we restrict this definition 
to require that there exist 2n — 1 globally defined functionally independent constants of 
the motion X{, % = 1, . . . , 2n — 1 that are quadratic in the canonical momenta pj. This 
clearly implies the relations 



^ V oxi dpi dpi dxi J 



2n- 1. 



The concepts of integrability and superintegrability also have their analogue in quan- 
tum mechanics. A superintegrable quantum mechanical system is described by 2n — 1 
(independent) quantum observables H — Xx, X%, . . . , X2 n -i that satisfy the commutation 
relations 

[H, Xi\ = HX - XiH = 0, i = 1, . . . , 2n - 1. 

The analogue of quadratic superintegrability in this case is that each of the quantum 
observables is a second order partial differential operator. Systematic studies of superin- 
tegrable systems have been conducted in spaces of constant curvature in two dimensions 

[3 1110013. 

In this article we solve the following problem. Given a Riemannian space in two 
dimensions with infinitesimal distance ds 2 = Ylf j = i9ij{u)du l du\ and u = {u l ,u 2 ), the 
classical Hamiltonian has the form 
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H = 9 l3 P l P, + V(u) 
and the corresponding Schrodinger equation is 

V 9 

where ^fg = det(<7y). Koenigs found all free Hamiltonians H = Y^Q^PiPj admitting at 
least two extra functionally independent constants of the motion of the form 



He obtained a number of families of solutions; in particular, spaces that admitted three 
extra quadratic constants. There must then be a functional relation between these and, 
furthermore, in each case there is a Killing vector, i.e., a function \i = Y^a=i al { u )Pi that 
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satisfies {H, fi} = 0. One of the three quadratic constants is a square of the Killing vector 
//. 

The problem we solve here is supplemental to that of Koenigs: Suppose we have a 
Hamiltonian H = Yl^PiPj + V( u ) that admits a Killing vector. We determine the 
potentials that correspond to superintegrability, i.e., potentials such that we can find at 
least two extra functionally independent quadratic constants of the form 

2 

A = ^2 " iJi -")l''l'.i + • 

A necessary condition that this be possible is that the Riemannian space be one of the 
four listed by Koenigs: 

1. ds 2 = (x + y)dxdy 



2. ds = - — + b dxdy 

\{x-y) 2 J 

3. ds 2 = (ae~ {x+y)/2 + be~ x ~ y ) dxdy 



a e 



(x~y)/2 + e (y-x)/2\ + fe 



4. ds = — r — dxdy. 

( e (x-y)/2 - e (y-x)/2) 2 

The first of these spaces, type one, or D±, has been treated in detail in an earlier 
paper [V . Here we treat the remaining three Darboux spaces in a similar and unified way. 
SectionsEJini andEJare devoted to the spaces D 2 , D 3 , and -D 4 , respectively. In each space 
we follow the same pattern. 

1. We first consider a classical free particle system and give the free Hamiltonian H , 
the Killing vector K and the two Killing tensors X\ and X 2 in a space with a 
conformally Euclidean metric (real or complex). We choose coordinates u and v 
in which the first order constant is K = p v , hence u is an ignorable variable, not 
appearing in the metric or in the Hamiltonian. 

2. We present an embedding of the two-dimensional Darboux space into a three- 
dimensional flat space. 

3. We present a polynomial relation between the four integrals of motion H, K, X\ 
and X 2 , and also the polynomial algebra generated by these integrals. 

4. We consider the quantum mechanics of a free particle in the corresponding Darboux 
space, i.e., write the corresponding Hamiltonian and integrals of motion as linear 
operators. We then establish that the relations between these operators are the 
same as those between the classical quantities. 

5. We use the fact that the Killing vector K generates a one-dimensional Lie transfor- 
mation group to classify all integrals of motion 

A = aX\ + bX 2 + cK 2 (1.1) 



4 



into conjugation classes. Each class gives rise to a coordinate system in which the 
Hamilton- Jacobi and Schrodinger equations allow the separation of variables. We 
construct these separable coordinate systems explicitly and solve the corresponding 
separated equations (classical and quantum). 

6. By construction, the free classical and quantum systems in Darboux spaces are 
all quadratically superintegrable: they have 3 functionally independent integrals of 
motion. We introduce potentials that do not destroy this superintegrability. Thus 
we present systematically all superintegrable classical and quantum systems of the 
form 

H = H + V(u,v) (1.2) 

where H$ is the free Hamiltonian in the space D 2 , -D3, or D4. To obtain this result 
we make use of the fact that to be quadratically superintegrable, a Hamiltonian in 
a Darboux space must allow the separation of variables in at least two coordinate 
systems. 

A separate section, Section |3J is devoted to the relation between superintegrable sys- 
tems in Darboux spaces and two-dimensional spaces of constant curvature. 

2 Darboux spaces of type two 

2.1 The free particle and separating coordinate systems 

If we allow rescaling of the variables x and y, as well as the Hamiltonian H then we can 
always take H to be of the form 

In the coordinates x = |(u + iu), y — \{v — iu) this Hamiltonian becomes 

u 2 (p 2 u + p 2 v ) 

H =r- 



u 2 + l 

Associated with the Hamiltonian are three integrals of the free motion 



r v Mgg - u 2 p 2 u ) {v 2 -u A )p 2 v + u 2 (l-v 2 )p 2 u 

K =p v , X 1 = 7-— \-2up u p v , X 2 = =— \-2uvp u p v . 

u 2 + 1 u + 1 



These three integrals satisfy the following polynomial algebra relations: 

{K,X 1 } = 2(K 2 -H ), {K,X 2 } = X 1 , {X U X 2 } = AKX 2 . (2.2) 
They are functionally dependent via the relation 

X\ - AK 2 X 2 + 4^X2 - AH 2 = . (2.3) 
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The corresponding problem in quantum mechanics can be obtained via the usual 
quantisation rules and symmetrisation: 

u 2 2v 
H = —-(d 2 u + d 2 ), K = d v , X l= M - u 2 d 2 u ) + 2ud u d v + d v , 

u z + 1 (ir + 1) 

X 2 = —5— ((v 2 - u 4 )d 2 + u 2 (l - v 2 )dl) + 2uvd u d v + ud u + vd v --. 
u z + 1 V / 4 

where the constant in the last expression is taken for convenience. The commutation 
relations are identical with those of the corresponding classical algebra: 

[k, x x \ = 2(k 2 - h ) [k, x 2 ] = x u [x u x 2 ] = 2{k, x 2 }. 

Here {K, X 2 } = ^(KX 2 + X 2 K). The operator relation (that exists in analogy with the 
functional relation in the classical case) is 

X\ - 2{K 2 , X 2 } + AH X 2 - 4H 2 -H + AK 2 = 0. 

The line element ds 2 = (du 2 + dv 2 )(u 2 + l)/u 2 can be realised as a two-dimensional 
surfaced embedded in three dimensions by 

X = 2&±±, Y-T=^±l, r + r=-g^±^±j^±I-^csinh M , 
u u 8u 8 

in which case, 

ds 2 = dX 2 + dY 2 - dT 2 = ^±(du 2 + dv 2 ) . 

u 

We wish to determine all the essentially different separable coordinate systems for 
the free classical or quantum particle. In order to do this we need to consider a general 
quadratic constant of the form A = aXi+bX 2 +cK 2 . Under the adjoint action of exp(aK), 
Xi and X 2 transform according to 

X x — X 1 + 2a{K 2 - H ), X 2 ->X 2 + aX 1 + a 2 {K 2 - H ). 

From these transformation formulae we see that if b ^ we can always take A in the form 
A = X 2 + (3K 2 . If b = then there are two representatives possible: X\ or K 2 . We have 
the following cases: 

X 2 + f3K 2 , X u K 2 . (2.4) 

We now demonstrate the explicit coordinates for each of these representatives using 
methods of our previous article[l . 

2.1.1 Coordinates associated with X 2 + (3K 2 

If we choose (3 = b 2 , b ^ suitable coordinates u, if are 

u = b cosh uj cos ip, v = b sinh uj sin ip , (2.5) 
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the standard form of elliptical coordinates in the plane. The classical Hamiltonian has 
the form 

H = pI+p% 

sec 2 ip — sech 2 c<j + b 2 (cosh 2 uj — cos 2 </?) 
The corresponding quadratic constant, expressed in these coordinates is 

2 2 (sec 2 (f + b 2 sin 2 <£>) p 2 + (sech 2 u; — b 2 sinh 2 o>) p 2 

X2 ~\~ b K = -. 2 — \ 7 o" \ ■ 

(sec 2 (p — sech uj) + b 2 (cosh u — cos 2 (p) 

The Hamilton- Jacobi equation is 



dS\ 2 /dS_ 
duj J \dip 



sec 2 if — sech 2 c<j + b 2 (cosh 2 to — cos 2 <p) 
with solutions of the form 



E. 



where /3i + /3 2 = —X/Eb 2 , fa fa — — 1/& 2 , $ = cos 2 = cosh 2 c<j. The corresponding 
Schrodinger equation 

(d 2 +d 2 ) * 

1 /.vl/ 



sec 2 cp — sech 2 o; + b 2 (cosh 2 uj — cos 2 tp) 
has solutions of the form 

* = Vcoi^cosh^ (ismhcu,~Eb \ Ps™ ^sin^, 



J = 1,2 



where S™ (z, «) and Ps™(t, k) are spheroidal functions jH] and i? — — 1 



m — 



•r 



2.1.2 Coordinates associated with X 2 

Here we use polar coordinates: 

u = rcos9, t> = rsin#. (2.6) 

The classical Hamiltonian has the form 

„ r 2 p 2 + p 2 



r 2 + sec 2 9 

and the corresponding quadratic constant is 

r 2 sec 2 Op 2 — 



Xo 



r 2 + sec 2 9 
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The Hamilton- Jacobi equation in these coordinates is 

\or ) \89 

r 2 + sec 2 9 

with solution 



E, 



^ / J^ r 2 _|_ \ 

S(r, 9) = V Er 2 + A — vA arctanhy 



- VX log (V\sin9 + VE- A COS z (7 

1 /= , / (-£/ + A) cos 2 9 — 2E 
H — v hi arccosh ' 



2 V (E ~ A) cos 2 

The corresponding Schrodinger equation is 



2 d 2 # 8 ( 1 <9^ 

r ~ r- + sin 



dr 2 89 \sin# 89 

r 2 + sec 2 9 

and has solutions of the form 



Em 



# = Vr sin C £+ i (y/^E rj P^{cos 9) , E = m 2 -^, 
where C v {z) is a Bessel function and P™(cos9) is an associated Legendre polynomial [H|. 

2.1.3 Coordinates associated with X\ 

A suitable choice of coordinates is 

u = £ri, v = ^e-V 2 )- (2-7) 
The classical Hamiltonian in these coordinates has the form 

The corresponding quadratic constant is 



x 1 



e+v 2 + ji + ^ 



The Hamilton- Jacobi equation has the form 
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which has solution 



y/E? + E-\? A / \e-2E 

o[^,T]) = — - — — arctanh 1 



£ 2 2VE \2y/Ey/E£ 4 + E- A£ 2 

+ 4E log (yE(2E£ 2 - A) + 2E^E£ A + E - A£ 2 ' 



jE^fTETW A / \r] 2 + 2E 

arctanh 1 



r/ 2 2v^E \2\/iV££ 4 + £ + A£ 2 ^ 

+ 4E log (v^Ef 2 + A) + 2Ey/E£* + E + 
The corresponding Schrodinger equation is 

a?* + 



e 2 + ^ 2 + ^ + i 



Typical solutions are 



* = — M x ^VEe)M- x ,,(^V 2 ) 



where M X)li (z) is a Whittaker function and .E = 4/i 2 — 4. 
2.1.4 Coordinates associated with K 2 

The representative K 2 has associated with it the coordinates u and v , in which the ignor- 
able variable has a fundamental role to play. The Hamiltonian and constant associated 
with this separation have already been given. The Hamilton- Jacobi equation has the form 

« 2 /7a?Y + m"i=tf. 



u 2 + 1 \ \ du J \dv J 
which has solution, with separation constant c, 



u 2 (E — c 2 ) + E 

S(u, v) = ^fu 2 (E -c 2 ) + E - VE arctanh W — ^ + cv . 

The corresponding Schrodinger equation has the form 

It +1 

Typical solutions are 

m = ^C V (Vm 2 - E u) e mv 

where E = v 2 — |. 

It is no surprise that the Hamiltonian is separable in elliptic, parabolic and polar 
coordinates, since, if we write the classical equation H = E in u, v coordinates we obtain 
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This equation is essentially the same form as a flat space superintegrable system with 
Cartesian coordinates u,v and potential a/u 2 , viz 

pI + pI + ^-e = o. 

It is known to be solvable via the separation of variables Ansatz in elliptic, Cartesian, 
polar and parabolic coordinates. This correspondence between flat space superintegrable 
systems and their curved analogues is essentially the way all the curved superintegrable 
systems can be obtained and is discussed in more detail in section El 

2.2 Superintegrability for Darboux spaces of type two. 

In this section we address the problem of superintegrability for the Hamiltonian 

H = ^4±f>. (2.8) 
u z + 1 

This is done in exactly the same manner as it was for the Darboux space of type 1 in 
a previous paper pQ. The free space Hamiltonian is given and we compute the possible 
potentials that correspond to superintegrability. There are four possibilities: 



[A] H = -J^ U+pl + a, (±u 2 + vA +a 2 v + ^ 2 



A basis for the additional constants of the motion is 

«i ( 9 w 2 + 4t> 2 \ a 2 ( o 4v 2 \ 2a 3 v 
R 1 = X 1 + -±v(u 2 + — — + -f U + ^ 



2\ u 2 + l J 2 V u 2 + l u 2 + l 



R 2 = K 2 + a\v 2 + a 2 v . 



These, along with R = {Ri, R 2 }, form a quadratic algebra 

that is determined by the identity 

R 2 = 16R 3 2 - AaiR\ - 32HR 2 2 - 8a 2 R ± R 2 

+ 8a 2 ifi?i + 1Q(H 2 + a x H - a l a 3 )R 2 + Aa 2 2 H - 4aja 3 . 

The classical equation of motion H — E = is 

2,2, / 1 2 , 2 \ , . a 3 ~ E n 

Pu + Pv + a i ( + v J + a 2 v -\ — E — 0. 

The basic form of this equation is a superintegrable system in flat space, but with rear- 
ranged constants, which is solvable via separation of variables in Cartesian and parabolic 
coordinates. 
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This accords with the fact that the leading part of a quadratic constant for this Hamil- 
tonian will be an element of the orbits represented by Xi and K 2 . So this Hamiltonian 
also separates in the 'parabolic' coordinates £, r\ ()2.7|) and in these coordinates takes the 
form 

V\ + V\ + \a x (e + r/ 6 ) + \a 2 (£ 4 - V 4 ) + «s + £ 



H 



e+v 2 + jz + j F 



Adding the same potential and coordinate functions to the quantum Hamiltonian H 
and its corresponding commuting operators X\ and K 2 , we obtain the operators 



H = H + -^—^ ( ai (\u 2 + v 2 ) +a 2 v + ^§ 



j 2 

■u 2 + 1 V " V 4 J * 



Ri = X i + — -v u H - — — + — - \u + 



2 V w 2 + 1 / 2V u 2 + lj u 2 + l 
R 2 = K 2 + a\v 2 + a 2 v . 

Ri and R 2 commute with H and along with R = [R\, R 2 ], obey the corresponding quantum 
quadratic algebra relations 

[R, Ri] = -24R 2 2 + 4a 2 i?i + 32HR 2 - 8H 2 - 8ai# + 6a x + 8a ia3 , 
[R, R 2 ] = -4ai-Ri - 4a 2J R 2 + 4a 2 # 



and the operator identity 

R 2 = l6R^ 2 i -4a 1 R 2 1 -32HR 2 2 -4:a 2 {R 1 ,R 2 } 

+ %a 2 HRi + 16H 2 R 2 + lQ ai HR 2 - 4ai(4a 3 - ll)R 2 
+ 4(a 2 2 + 8ai) # - Ab 2 2 (a 3 + |) . 

[B] H=^L 2 u+P 2 v + b 1 (u 2 + v 2 ) + h ^ + h 4). 

The additional constants of the motion have the form 

Ri = X 2 + Hl±j£ (fc ( M 2 + v 2 ) -b 2 -b 3 ^), R 2 = K 2 + M 2 + - 2 • 

The corresponding quadratic algebra relations can be determined, using (j2.9J) . from the 
identity 

R 2 = l6R 1 R 2 2 -l6b 1 R 2 -l6HR 1 R 2 + 32b 1 (H-b 2 -b 3 )R 1 + l6(H + b 3 -b 2 )HR 2 
- 16(6! + b 3 )H 2 + 32^(62 - b 3 )H - I6&1O2 - &s) 2 • 

The equation of motion H — E = becomes 
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This is a superintegrable system in flat space, but with rearranged constants, which is 
solvable via separation of variables in Cartesian, polar and elliptic coordinates. Again, 
this agrees with the observation that for this Hamiltonian we have quadratic constants 
with leading parts K 2 , X 2 and X 2 + (3K 2 . In the latter two coordinate systems, the 
Hamiltonian takes the forms: 

(i) Elliptical coordinates ()2.5j) 

vt + V% + i6i6 2 (sinh 2 2uj + sin 2 2ip) + 6 2 (sec 2 ip — sech 2 cj) + 6 3 (cosec 2 y? + cosech 2 a;) 



H 



sec 2 if — sech 2 a; + b 2 (cosh 2 to — cos 2 ip) 
(ii) Polar coordinates ()2.6j) 

r 2 p 2 + + hr A + b 2 sec 2 9 + 6 3 cosec 2 6 l 
r 2 + sec 2 9 

The corresponding quantum algebra relations are 

= -8{£i, 4} + 8Mi + 124- 8# 2 + 8(6 2 - b 3 - \)H, 
[R,R 2 ] = 8R 2 2 - I661E1 -8^2 + 1661^-1661(63 + 63 + |), 
i? 2 = 8{i?i, i?^} - 8fl"{^i, R 2 } + 16# 2 i? 2 - 166ii? 2 - 76R 2 2 

+ 326iMi - 861(4(63 + 6 2 ) + 3)i?i + 16(6 3 - 6 2 + f)HR 2 

- 16(6i + 63 + \)H 2 - 861(4(63 - 6 2 ) + 3)H + 61 (36 + 486 3 - (4(6 3 - 6 2 ) + 3) 1 

[C]H= e + v> + ^ + ^ • 

The additional constants of the motion are 

Cie 2 (r/ 4 + l) + C2(r ? 4 + l)-c 3 (e 4 + l) 



Ri = X x + 
R 2 = X 2 + 



(eV + l)(e 2 + r/ 2 ) 
Ci(e 2 + ^)-C2(r/ 4 -l)-c 3 (^ 4 -i; 



4(fV + 1) 

The corresponding Poisson algebra can be determined from the identity 

R 2 = AR\R 2 - (c 2 + c 3 )Rl + 16HR 2 2 - A Cl RiR 2 + 2ciC 3 i?i - 16H 2 R 2 
+ 4(c 2 + c 3 )H 2 + (c 2 - Ac 2 c 3 )H - c 2 c 3 

The Hamiltonian can be written in separable form for the following coordinate systems. 

(i) Displaced Elliptic coordinates £ = 6' cosh u' cos ip', 77 = 6' sinh u/ sin <£>'. 

p 2 , + p 2 , + Ci6 /2 (cosh 2 u/ — cos 2 if') + c 2 (sec 2 <// — sech 2 c</) + c 3 (cosec 2 <// + cosech 2 ci/) 
6 /4 (cosh 4 a;' — cos 4 ip 1 — cosh 2 a;' + cos 2 if') + sec 2 <// + cosec 2 ^' + cosech 2 u/ — sech 2 u/ 
These coordinates are not those given in ()2.5|) and are related to u and v by 

u = -b' 2 sinh 2u/ sin 2</?' , v = -b' 2 (cosh 2a;' cos 2<f' + 1) . 
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(ii) Polar coordinates £ = r' cos 8', 77 = r'sin9'. 

r /2 p 2 , + p 2 ,, + cir' 2 + c 2 cosec 2 6' / + c 3 sec 2 6' 
r' 4 + sec 2 9' + cosec 2 #' 

These coordinates are not those given in (j2.fi |) and are related to u and i> by 

m = -r' 2 sin 20' , u = -r' 2 cos 20' . 
2 ' 2 

The corresponding quantum algebra relations are 

[R, Ri] = -2R\ - 2 Cl R! - I6HR2 + 8H 2 - 6H 
[R,R 2 ] = 2{R u R 2 }-(c 2 + c 3 )R 1 -2c 1 R 2 + c 1 c 3 

R 2 = 2{Rl R 2 } + 16HR 2 2 - (c 2 + c 3 + 4)R 2 -2c 1 {R 1 , R 2 } + 2c 1 (c 3 + 2)R 1 

- 16H 2 R 2 + \2RR 2 + 4(c 2 + c 3 )H 2 + (c 2 - 4c 2 c 3 - 3(c 2 + c 3 ))# - ±(3 + 4c 3 )c 2 

The equation of motion if — = is 

This is a superintegrable system in flat space, but with rearranged constants, which is 
solvable via separation of variables in Cartesian, polar and elliptic coordinates. 

rril rr U 2 {p 2 u +p 2 v + d) 

The additional constants of the motion are 



2dv d(u 2 + v 2 ) 

R l = X l + —— ) R 2 =X 2 + 1 \ \ K = Pv . 
u l + 1 u z + 1 



The corresponding Poisson algebra relations are 

{K,R 1 } = 2K 2 -2H + 2d, {K,R 2 } = R 1 , {R 1} R 2 } = -AKR 2 . 

The functional relation between these constants is 

R\ - AK 2 R 2 + A(H - d)R 2 - AH 2 + AdH = . 

The Hamiltonian can be written in separable form for all the possible types of separable 
coordinates we have discussed viz. 



(i) Elliptic coordinates (|2.5|) . 



H 



Pt + vt, + & 2 c?(cosh 2 uj — cos 2 (f) 



6 2 (cosh uj — cos 2 ip) + sec 2 tp — sech u> 
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(ii) Polar coordinates ([2.6)1 . 

_ r 2 p 2 + + <ir 2 
r 2 + sec 2 9 

(iii) Parabolic coordinates (|2.7|l . 

p|+p 2 + rf(e 2 + ^ 2 ) 



The corresponding quantum algebra relations have the form 

[AT, fli] = 2K 2 - 2H + 2d, [K, R 2 ) — R\, [R u R 2 ] = 2{K, R 2 }. 
The operator identity satisfied by the defining operators of the quantum algebra is 

R\ - 2{K 2 , R 2 } + AHR 2 - AdR 2 + AK 2 - AH 2 + (Ad - 1)H = . 
The equation of motion H — E = is 

p 2 u +p 2 v + d-E-^ = 0. 

u 

This is a superintegrable system in flat space, but with rearranged constants, which is 
solvable via separation of variables in Cartesian, polar, elliptic and parabolic coordinates. 

3 Darboux spaces of type three. 

3.1 The free particle and separating coordinate systems 

With rescaling and translation of the variables x and y the Hamiltonian H has the form 

H ° = 1 + e -(^)/ 2 PxPy - (3 ' 1} 

In coordinates x = u — iv, y = u + iv we can write this Hamiltonian in positive definite 
form 

l e 2u (p 2 u +p 2 v ) 
A e u + l 

Associated with the Hamiltonian are three integrals of the free motion 

1 e 2u 2 le"(e u + 2) 2 1 u . 

K = p v , X 1 = - — - - cosvp u - - ^ cosvp v + -e smvp u p v , 

1 e 2u . 2 le u (e" + 2) . 2 1 n 
2 = 4^TT Sm?;Pw ~ 4 e" + l Smt; ^ ~ 2 6 C0S ^"^ • 
The integrals satisfy the polynomial algebra relations 

{K, Xx} = -X 2 , {K, X 2 } = X h {X h X 2 } = KHq. 
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They are functionally dependent via the relation 

X\ + X\ - Hi - H K 2 = 0. 

The corresponding problem in quantum mechanics can readily be obtained via the 
usual quantisation rules and symmetrisation. 

1 e 2u 

H ° = Z^iW + %)> K = d v , 



1 e 2u le u (e u + 2) 1 1 1 

Xi — cos vd 2 — - cos v d 2 H — e u sin vd u d v + -e u cos vd u + -e u sin vd v , 

4e u + 1 4 e u + 1 2 4 4 

. 1 e 2 « 



sinvd 2 — — — - — -t— sint><9 2 — -e" cosvd u d v + -e u sint><9„ — -e u cosvd v 
4 e u + 1 4e" + l 2 4 4 



The commutator algebra obtained has the same form as the Poisson algebra, and the 
identity relating the operators is 

X\ + X 2 - Hi - H K 2 + ^H = 0. 

The line element ds 2 = (e~ u + e~ 2u )(du 2 + dv 2 ) can be realised as a two-dimensional 
surface embedded in three dimensions by 

X = vVe~ u + e~ 2u , Y -T = Ve~ u + e^ 2 " 

Y + T = (l-v 2 )Ve~ u + e- 2u + log (l + 2e _ " + 2\fe~ u + e~ 2u j + ^arctan ^2 \/e~ u + e- 2u j , 
in which case, 

rfs 2 = rfX 2 + rfF 2 - dT 2 = (e' u + e- 2u )(du 2 + dv 2 ) . 

Just as we have done in other cases, we wish to determine all the essentially different 
separable coordinate systems for the free classical or quantum particle. To do this we 
need to consider a general quadratic constant of the form A = aX\ + bX 2 + cK 2 . Under 
the adjoint action of exp(aK), X\ and X 2 transform according to 

X\ — > cos a X\ — sin a: X 2 , X 2 — > sin a; X\ + cos a; X 2 . 

From this transformation law we see that A can take five different forms: 

K 2 , X u X 1 + 1 K 2 , X! + ?X 2 , X ± + iX 2 -K 2 . (3.2) 

We now demonstrate the explicit coordinates in the case of each of these representa- 
tives. 
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3.1.1 Coordinates associated with K 2 

These are the coordinates associated with the ignorable coordinate v and the Hamiltonian 
has already been given in the u, v coordinates. The Hamilton- Jacobi equation is 

4e u +l \[du) h \~dv) ~ L 

with solutions 



n , , y/4E(l + e u ) - c 2 e 2u y/E , / ^E(e u + 2) 
b(u,v) = — - arctanh 



ce u c V JAE(1 + e u ) - c 2 e 2u 



+ i log (i(c 2 e u - 2E) + cy/AE(l + e u ) 



c 2 e 2u l + cv . 



The corresponding Schrodinger equation is 

1 e 2u 



Ae u + 1 

with solutions of the form 



(d 2 + d 2 )^ = B9, 



* = e-«/ 2 M_ 1/v ^ i±m (4 v / =^e-«)e 
3.1.2 Coordinates associated with X\ 

For the second representative in (|3.2|) . a suitable choice of variables is 

f = 2e- u/2 cos | , 77 = 2e" M/2 sin ^ . (3.3) 
In terms of these coordinates the classical Hamiltonian has the form 

-"0 " 



4 + e + v 2 

and the corresponding quadratic constant is 

(2 + r/ 2 )pf - (2 + e)p 2 



n 



2(4+e+^ 2 ) 

In £, 7] coordinates the classical Hamilton- Jacobi equation is 



A + e + V 2 

which has the solution 



E. 



S = \^Ee + 2E-X+ (^j=f) log (E + y/E? + 2E-\) 
+ ^vV E V 2 + 2E + X + J log (E + V E V 2 + 2E + 



A 
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The Schrodinger equation is 

(gg+g); 

4 + e + r] 2 

which has typical solutions 



Em 



* = ^(A- 2 i 3 )/v^( ± ( 4 ^) 1/4 ^-(A + i ; )/v^( ± ( 4 ^) 1/ ^ 

in terms of parabolic cylinder functions D v (z) [HUH]. 
3.1.3 Coordinates associated with X\ + jK 2 

For the third case it is convenient to take the representative as b 2 Xi + 2K 2 . Here we 
identify coordinates via 

£ = b cosh uj cos ip , r/ = 6 sinh uj sin </? . (3.4) 
The classical Hamiltonian has the form 



pi+pl 



26 2 (cosh 2u - cos 2ip) + |6 4 (cosh 2 2u - cos 2 2ip) 
and the corresponding quadratic constant in these coordinates is 

(8cos2<^ — b 2 sin2ip)p 2 + (8cosh2co> + fe 2 sinh2u;)p 2 

b Xi + 2i^ = 



86 2 (cosh 2u - cos 2^) + 6 4 (cosh 2 2uo - cos 2 2<p) 
In the a; coordinates the classical Hamilton- Jacobi equation has the form 



dS\ 2 fdS_V 
du J \dtp J 



26 2 (cosh 2u - cos 2ip) + ±6 4 (cosh 2 2u - cos 2 2ip) 
and has the solution 



E 



where ai + a 2 = — P\ — /3 2 = —8/b 2 , ai« 2 = — /3i/? 2 — 4:X/Eb 2 , Q = cosh2u;, $ = cos2<^. 
The corresponding Schrodinger equation 

<9 2 # + d 2 ^ 



26 2 (cosh2cu - cos 2<p) + ±6 4 (cosh 2 2u - cos 2 2<p) 
separates with \1> = §((p)Q(u) in the equations 

(<9j + 2& 2 £cos2<^ + §6 4 £cos4<^ + A + §& 4 £)$ = 
(dl - 2b 2 E cos2u -\tfE cos - \-\b A E}VL = 
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which has typical solutions 



^2 



i 




with corresponding separation constant given by X m = fi m b 2 (l + b 2 )E/8. The functions 
appearing here are even and odd Whittaker Hill functions [TU] . 

3.1.4 Coordinates associated with Xi + iX 2 

In the case of a system specified by the fourth representative there are, in fact, no separable 
coordinates. However, in the coordinates 



the classical Hamiltonian takes the form 



H 



tyxPy 



2y-V3 + x 



and the corresponding constant is 



X, + iX 2 = 2yH - p 2 x . 



The solution of the Hamilton- Jacobi equation 



as as 

dx dy 



E 



2y~V 2 + x 



is 




The corresponding Schrodinger equation is 



2dM 



E^ 



2y- 1 / 2 + x 



which has solutions 
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3.1.5 Coordinates associated with X 1 + iX 2 — K 2 

In the case of a system specified by the fifth representative an appropriate choice of 
coordinates is 

2JJw \2JjIu J 



The corresponding classical Hamiltonian has the form 

h 2 pI - v 2 p 2 v 



(ji + v)(2 + fi-v)' 
and the quadratic constant is 



2 _ u 2 (fi + 2) fj.pl -fi 2 (u- 2)vpl 



X, - iX 2 - K 

{(J, + v) (2 + ix - v) 

In the fi, v coordinate system the classical Hamilton- Jacobi equation has the form 



,/SSV 2 f9S\ 2 



+ I/) (2 + /i - V 

which has the solution 

S(jjl,u) = ^/Efi 2 + 2Efi + A + \[E log (y/E{l + li) + v 7 ^ 2 + 2£/i + A 

- arctanh f , A + ^ 

\v / A v / ^A t2 + 2 ^/ i + A , 

+ VEv 2 - 2Ev + A + v 7 ^ log ^-/E(l - v) + VEv 2 - 2Ev + A 

— vA arctanh 



\ — Ev 



VX^Eu 2 - 2Ev + A 



The Schrodinger equation 



(/i + z/)(2 + /i-z/) 
separates with \1/ = A(/j,)B(u) into the equations 

(udnbidj - Ell 2 - 2£/i - p 2 ) = , (vd v {vd v ) - Ev 2 - 2£z/ - p 2 ^B(u) = 

and has solutions, in terms of the Whittaker function M\ tX , of the form [§] 

l=M^ p (2VE h) M_^ p (2VE 



'LLV 

If we write the classical equation H = E in £, rj coordinates we obtain 

pl+p 2 v -E(A + e + V 2 ) = 0. 

This is in the form of a flat space superintegrable system which can be solved by separation 
of variables in Cartesian, polar, hyperbolic and elliptic coordinates. 
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3.2 Superintegrability for Darboux spaces of type three. 

In this section we address the problem of superintegrability for the Hamiltonian 

H=- {Pu+p v> . (3.6) 
4 e u + 1 v ' 

We arrive at three possibilities: [A], [B], [C]. 

The additional constants have the form 

2^(2 + r/ 2 ) - 2a^(2 + £ 2 ) + a 3 ( V 2 - £ 2 ) 



i?i = X 1 + 

i?2 = -^2 + 



4(4 + e + 

a^^a - £ 2 + 4) + a 2 £(e - V 2 + 4) - 2a 3 ^ 



4(4 + e 2 + V 2 ) 

The corresponding quadratic algebra can be determined from the identity 

R 2 = HR\ + HR 2 2 + \(a 2 2 - a\)R l -\a l a 2 R 2 

- H 3 + \a z H 2 + ^(2a 2 + 2a\ - a\)H - ^a 3 (a 2 + a 2 ) . 

This Hamiltonian separates in a family of coordinate systems obtained by translating the 
given separable system via £ — > £ + a, 77 — > 77 — a. The corresponding quantum algebra 
relations are 

[i?, = —HR 2 + |cna2 , -R2] = -f^-Ri + Yg( a 2 ~~ a i) > 

-R 2 = HR 2 + HR 2 + |(o 2 — a 2 )i?i — ^a,\<x 2 R 2 

- # 3 + |(a 3 + |)# 2 + ^(2a 2 + 2a 2 , - a 2 ,)// - ^a 3 (a 2 + a 2 ) 

As in the case of free motion, the equation H = E becomes 

p\ + V% + aif + + a 3 - ,5(4 + £ 2 + r/ 2 ) = . 
Again, this is a superintegrable system in flat space but with rearranged constants. 

[B] H = — 4+e + v 2 — ' 

The additional constants are 

Ri = Xi + , (3 . 7) 
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[R, Hi] = 2R\ - AHR 2 - 2H 2 + (6 3 + §) # - |^ 



The corresponding quadratic algebra relations are determined by 

-R 2 = -AR\R 2 -(b 1 + b 2 )Rl + 4HR 2 2 + 2(b 1 - b 2 )HR 1 + \h{b 2 -b 1 )R 1 
+ AH 2 R 2 - 2b 3 HR 2 + \b\R 2 - {b x + b 2 )H 2 
+ {Ihih + b 2 ) - hb 2 ) H - ±bl(b! + b 2 ). 

This Hamiltonian separates in all the separable coordinate systems given in Section 12.11 
The Hamiltonian has the explicit forms 

(i) In u, v coordinates, 

_ e 2u (p 2 u +p 2 v + \bx sec 2 f + ±6 2 cosec 2 f + b 3 e~ u ) 
' ~ A{e u + 1) 

(ii) In the elliptical coordinates (|3.4j) . 

vt + V% + 6i (sec 2 (y? — sech 2 co>) + 6 2 (cosec 2 ^ + cosech 2 a;) + 6 3 6 2 (cosh 2 u — cos 2 <£>) 
26 2 (cosh 2uo - cos 2y>) + ±6 4 (cosh 2 2u - cos 2 2^) 

The corresponding quantum algebra relations have the form 

±h 2 

8°3 

[R,R 2 ] = -2{Hi, J R 2 }-(&i + fo 2 + l) J Ri + (6 1 - 6 2 )# + f(6 2 -6i)6 3 , 
-R 2 = -2{i? 2 ,i? 2 }-(6 1 + 6 2 + 5) J R 2 + 4M 2 

+ 2(6! - b 2 )HR! + 63(62 - 6i)Hi + AH 2 R 2 - (26 3 - l)HR 2 + \b 2 3 R 2 

- (61 + 6 2 - 2)H 2 + (|(6 3 + §)(&! + 6 2 ) - 63 - 6 X 6 2 -\)H- ^6 2 (6 X + 6 2 - 2). 

As in the case of free motion, we observe that equation H = E becomes 

pl + P 2 v + ^ + ^ + b 3 -E(A + e + V 2 ) = 0- 

This is a sup erintegr able system in flat space, with rearranged constants, that separates 
variables in Cartesian, polar and elliptic coordinates. 

rc , R ^-^ + Cl(/i + ^)+C2^+C 3 ^ 

[ J {n + u){2 + ii-u) 

The additional constants of the motion have the form 

Cj/iV + C 2 /iZ/ + 2C 3 (1 + /X - 2 //-I/ (/J - Z/) 2 

ill = Ax + «A 2 — — K 2 = K - C 2 C 3 — — . 

/iZ/(2 + /J, — V) jjLU jJL z V z 

The corresponding quadratic Poisson algebra relations can be determined from 
R 2 = -AR 2 R\ + Sc 2 HR l -Ac l c 2 R l + lQc 3 HR 2 



+ 16c 3 H 2 + 4(c 2 - 4cic 3 )# + 4c 2 c 3 . 
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The quantum algebra relations are 

[R, Ri] = 2R\ - 8c 3 H , [R, R 2 ] = -2{R U R 2 } - R x + Ac 2 H - 2c x c 2 , 

R 2 = -2{R 2 , R 2 } + &c 2 HR x + l6c 3 HR 2 - hk\ - ^c x c 2 R x 
+ 16c 3 H 2 + 4(c 3 + c 2 2 - Ac x c 3 )H + 4c 2 c 3 . 
As in the case of free motion, equation H = E becomes 

a superintegrable system in flat space with rearranged constants, that separates variables 
in polar and hyperbolic coordinates. 

_ /j, 2 pl - u 2 pl + d 1 fi + d 2 u + d 3 (/j, 2 + v 2 ) 
[ 1 ~ {ii + v){2 + n-v) 

The additional constants of the motion have the form 

liv(d x {y - 2) + d 2 (fi + 2) + 2d 3 (y - /i + /jlu)J 



Ri = X 1 + iX 2 -K 
R 2 = X 1 -iX 2 - 



+ z/)(2 + fi-v) 
(fji - v) [{ji - v){d\\i + <kv) - 2d 3 (fx 2 + v 2 + fjiv(2 + /x - v)) 
AfivQi + v)(2 + n - u) 
The corresponding quadratic Poisson algebra can be determined from 

R 2 = 4RiR 2 2 - AHRtRz + d\R x - AH 2 R 2 + 2(d 1 + d 2 )HR 2 - d x d 2 R 2 

+ 4# 3 - 2(di + d 2 )H 2 + \((di + d 2 ) 2 + d 3 (d 2 - difjH- d 3 {d\ - d 2 ). 

This classical system also separates in elliptical coordinates obtained by choosing new 
variables defined by the roots of the characteristic equation of Ri + R 2 , that is, the 
elliptical coordinates (|3.4|) with b = 2i. In these variables the Hamiltonian has the form 

H Pl+Pl + 2(c?i + d 2 )(cos2(f - cosh2cu) + 2(di - d 2 )(2ism2(p + sinh2cu) + 2d 3 (sinh4u; + 2isin40) 

8(cos 2tp — cosh 2uS) + 4(cosh 2 2uj — cos 2 2(p) 

The corresponding quantum algebra relations are 

= -2{R 1 ,R 2 } + 2HR 1 + R 2 + 2H 2 - [d 1 + d 2 + \) H + \d^ 2 
[R, R 2 ] = 2R\ ~~ 2HR 2 + \d\ 

R 2 = 2{R U R 2 2 } - hh\ - 2H{R U R 2 } + d 2 ^ - AH 2 R 2 + (2di + 2d 2 + h)HR 2 - d!d 2 R 2 
+ 4H 3 - (2d 1 + 2d 2 + l)H 2 + (i(d! + d 2 ) 2 + d 3 (d 2 - dt)) H - \d 3 (d 3 -d\ + d\) 

As in the case of free motion we observe that equation H = E becomes 
Pl + pl + d 1 + d 2 -4d 3 ■ (*" d 0(e-^) , Sda^ + ir,) 



v / ^-^) 2 + 4 >/e =^+4 (e - iv + >/g - ^ + 4) : 
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= (E-d 3 )(A + e + V 2 ), 

a superintegrable system in flat space with rearranged constants that separates variables 
in elliptic and hyperbolic coordinates. 

The additional constants of the motion are 

c ri 2 — £ 2 c £w 

D Y _L ' S D Y s ' 

-"■1 ~~ A l + T . , AO , o ; ti-1 — -^2 — ~ 



44 + e + ?7 2 ' 24 + ^ 2 + ?7 2 ' 

and X. The corresponding Poisson algebra relations have the form 

{K,R 1 } = -R 2 , {K,R 2 } = R U {R 1 ,R 2 } = HK, 
and the functional relation between these constants is 

2 

R\ + R 2 2 - HK 2 — H 2 + -H — — = . 
11 2 16 

This Hamiltonian separates in all of the four types of separable coordinate systems 
available, and the corresponding expressions for the Hamiltonian can be deduced from [2] 
by taking b 3 = c, b\ = b 2 = 0. 

The quantum algebra relations are 

[K,R 1 ] = -R 2 , [K,R 2 ] = Ri, [Ri, R2] — HK , 
and the associated operator identity is 

4 Darboux spaces of type four 

4.1 The free particle and separating coordinate systems 

With rescaling of the variables x and y, the Hamiltonian H can be taken in the form 

(e x ~ y — e y ~ x ) 2 

H ° = x-vi v-x! P-Py ( AA ) 

e x y + e y x + a 

In coordinates x = v + iu, y = v — iu, we can write the Hamiltonian as 

sin 2 2u(p 2 u +p 2 v ) 

riQ = . 

2 cos 2u + a 

It admits constants of the motion 

K = p v Xi — e 2v (—H + cos 2u p 2 u + sin 2u p u p v ) , X 2 = e~ 2v (—H + cos 2u pi — sin 2u 
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These integrals satisfy the polynomial algebra relations 

{X, X 1 } = 2X h {X, X 2 } = -2X 2 , {X 1; X 2 } = -8X 3 - AaKH . 
They are functionally dependent via the relation 

X X X 2 - X 4 - aK 2 H - Hi = 0. 
The corresponding quantum operators are 

A — sin 2 2u / a \ 

H ° = s 7T— i d l + d D , X i = e2v -Ho + cos 2u(d 2 v + d v ) + sin 2u(d u d v + d u ) ) , 

2 cos 2u + a V / 

X = d v , X 2 = e~ 2v (-H + cos 2u{d 2 v - d v ) - sin 2u(d u d v - . 

Their algebra is determined by the relations 

[K, Xi] = 2X 1? [X, X 2 ] = -2X 2 , [Xi, X 2 ] = -8X 3 - 4aKH - 4X, 

and the operator identity is 

-{X x , X 2 } - X 4 - aH K 2 - 5X 2 - Hi - aH = . 

The line element <is 2 = (2cosw + a)(du 2 + <iv 2 )/ sin 2 2m can be realized as a two- 
dimensional surface embedded in E(2, 1) by (assuming a > 2) 

X = \J a + 2 cos 2m>, V — T = Vet + 2 cos 2m, 

(a -2) . . /a - 2 2 . , /a - 2 2 , „ 

where 

/(a + 2) (cos 2m + 1) 2 
y 2(a + 2cos2w) V a + 2 

and II is an elliptic integral of the third kind [Hj. Then ds 2 = dX 2 + dY 2 — dT 2 . 

Just as we have done in other cases, we wish to determine all the essentially different 
separable coordinate systems for the free classical or quantum particle. To do this we 
need to consider a general quadratic constant of the form A = aX\ + 6X 2 + cK 2 . Under 
the adjoint action of exp(aX), X\ and X 2 transform according to 

Xi -> exp(-2o;)X 1 , X 2 -> exp(2a)X 2 . 

If we regard two such quadratic expressions as equivalent if they are related by a combi- 
nation of group motions and the discrete transformation observed above, then the equiv- 
alence classes of these expressions can be chosen to have the following representatives: 

X 2 , X 2 , 7 X 2 + X 2 , X! + X 2 + 7 X 2 . (4.2) 

In the last of these are three cases to distinguish: 7 = 0, 7 = 2 and 7 7^ 0,2. The various 
separable systems involved can now be computed. 
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4.1.1 Coordinates associated with K 2 

These are the coordinates associated with the ignorable coordinate v and the Hamiltonian 
has already been given in the u, v coordinates. The Hamilton- Jacobi equation is 

^2u ((BS\> B . 



2 cos 2m + a \ \ dv J \du 

It has typical solutions 



S(u,v) = — i\ogyi(c cos 2m — E) + cy E(a + 2 cos 2u) + c 2 sin 2m 

1 ^ , / E(a + 1) + c 2 + (E - c 2 ) cos2w 

+ — v 7 ^ a + 2 arctanh v ; v ; 

2c \ v/£(a + 2) (£(a + 2 cos 2m) + c 2 sin 2 2m) 



1 ^ r , / E(a- 1) + c 2 + (£ + c 2 )cos2m 

+ — a/£ a- 2 arctanh v y v ' =. 

2c V ^(a - 2) (£(a + 2 cos 2m) + c 2 sin 2m) 



The corresponding Schrodinger equation is 

sin 2 2m fd 2 ^ d 2 ^\ 
1 = E^ , 

2 cos 2m + a \ dv 2 du 2 J 

which has the solution 

* = 2^1 Q(A - e+ - e_), ^(A + e+ + e_), e+ + i sin 2 u ) e Xv , 

where 



and 2-F1 is a Gaussian hypergeometric function [TT] . 

4.1.2 Coordinates associated with X 2 

If we choose new coordinates 



2 vr , j , °\2 
then the Hamiltonian takes the rational form 



4uV(p 2 +p 2N 

-"0 " 



6± = - + -v / l-(«±2)^, (4.3) 



x = lo S ( - ^) J > 1/ = log f^(/i + ^) ) , (4.4) 



(a + 2)u 2 + (a-2)z/ 2 

In this case the corresponding choice of coordinates has already been given, and the 
quadratic constant in these coordinates is 

4(« + 2)u 2 p 2 -4(a-2)z/ 2 p 2 
2 (a + 2)u 2 + (a - 2)z/ 2 ' 
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The Hamilton-Jacobi equation 

1/<V I I — I + I — 



as\ 2 /as x2 ' 



9/i / \ <9i/ 



(a + 2)/i 2 + (a - 2)z/ 2 

has solution 



S(fi,v) = i^{a - 2)E + Xfi 2 - iy/(a - 2)E arctanhW ^ 2 ^ ^ A/X 



+ v/(a + 2)E - Az/ 2 - v / (a + 2)£ arctanh, / ^ y 2 ^ oA ^ 

y (a + 2j£/ 

The corresponding Schrodinger equation has Bessel function solutions of the form 

4.1.3 Coordinates associated with 7X2 + K 2 

In the case of the third representative the transformation 

/i = c cosh uj cos ip , z/ = c sinh sin y9 (4.5) 
gives the classical Hamiltonian 

H= KpI + pD 

(a — 2)(sech 2 a; — sec 2 </?) — (a + 2) (cosech 2 c<j + cosec 2 (/?) 
The classical constant associated with this coordinate system is 
r 2 

4 2 

((a — 2) sec 2 </? + (a + 2)cosec 2 y?) p 2 + ((a — 2)sech 2 a; — (a + 2)cosech 2 u;) p 2 
(a — 2)(sech 2 a; — sec 2 <p) — (a + 2)(cosech 2 u; + cosec 2 y?) 

The Hamilton-Jacobi equation in these coordinates is 

4(P 2 +P 2 ) 

2 2 — 

(a — 2)(sech — sec 2 if) — (a + 2)(cosech a; + cosec 2 v?) 

and has solutions of the form 

S(u,<p) = i log (^/\ (A cos 2f + 2aE) + \^8E + 4aE cos 2y? - A sin 2 2y^j 

1 /7 ^ 1 ( 2Ea + 8E- A + cos2y?(A + 2£a) \ 

- -v/ a + 2 £ arctanh — yv ; = 

4 \ 2 V 7 ^ + 2)E(8E + AaE cos 2y? - A sin 2 2(p) J 
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1 n ^ i ( -2Ea + 8E - X + cos2(p(-X + 2Ea) \ 

- -J(a - 2)E arctanh — ^ v ; — 

4 \2y/(a-2)E(8E + 4aE cos2y? - \sm 2 2ip) J 

+ - VX log (y/\ (A cosh 2u - 4E) + A ^A sinh 2 2cj - 8£ cosh 2u - iaE^j 

, / AE + A + AaE + cosh 2ujUE - A) 

+ v a + 2 arctan ( — ^^^^^^^^^^^=^^^^^^^^^^^= 

2 J (a + 2)E(X sinh 2 u - 8E cosh 2cu - AaE) 
-4E + A + 4a£ + cosh 2u{AE + A) 



+ \J a — 2 arctan 



2\j[a-2)E[X sinh 2 - 8£ cosh 2cj - 4a£) 
The Schrodinger equation has the form 



,'d 2 ^ a 2 * 

4 ( 1 

du 2 dip 2 



E^ 



(a — 2) (sech 2 a; — sec 2 <p) — (a + 2) (cosech 2 co> + cosec 2 <^) 
with corresponding solutions 

\& = (sin ip sinh u;) e ~ (cosy? cosh oj) t+ 2-F1 ^ + , + - — , e_ + -, sin 2 cp 

_, / e + + e_ - A e + + e_ + A 1 . 2 

x 2F1 ( 2 ' 2 ' C 2' _ ^ 

with e± defined by (|4.3|) . 

4.1.4 Coordinates associated with Xi + X 2 + 7-f^ 2 

For the coordinates corresponding to the fourth representative we make the transformation 
u = arctan(exp a), v = (3/2, so our Hamiltonian has the form 

g= _ 4 ^ ±f eAV| 
a — 2 tanh a 

This can be realised in terms of projective coordinates on a two-dimensional complex 
sphere via si = cosh a cosh /3, s 2 = z cosh a sinh /3, s 3 = i sinh a where s 2 + s 2 + s 2 = 1. 
The Hamiltonian can be written as 



H = A 



^1 ~j~ ^2 ~j~ 



These two ways of realising the classical Hamiltonian are useful in determining the various 
possible separable coordinate systems. 

We consider the most general case first, i.e., 7 ^ 0,2. We make use of the transfor- 
mation equations 

■ XY + 1 

sinh a = 1 . , 

2\[XY 
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tanh j3 



2^{A + X + A_){A_X + A + )(A + Y + A„)(A_Y + A, 



(A\ + A 2 _ )(XY + 1) + 2A + A^(X + Y) 
applied to (|4.6|) to give classical Hamiltonian in the form 

.X(A + X + A„){A^X + A + )p 2 x - Y{A + Y + A_){A_Y + A + )p 2 Y 



H = -16XY- 



A + A_ {X - Y) {[a + 2)XY - a + 2 



The corresponding classical constant associated with this coordinate system is 
Al + A 2 „ 

^ + X 2 + 



16 



Vx 



(A + X + + A+) (a(A+Y + A_)(A_Y + A+) - 2A + A_(Y 2 - 1)^X 

A + A_ {A% - A 2 _)(X - Y) ((a + 2)XY - a + 2) 
(A+F + A_){A_Y + A + ) (a(A + X + A_)(A_X + A + ) - 2A + A_{X 2 - 1))y 

1 R V / 



Vy 



A + A_ {A\ — A 2 _){X — Y) ({a + 2)XY -a + 2 
The Hamilton- Jacobi equation has the form 

X(A + X + A^)(A^X + A+) (|Q 2 + Y(A + Y + A_)(A_Y + A+) ' 
A + A„ (X - Y) {{a + 2)XY -a + 2 

and solutions 



K ] 4A^A1\ X } X\\{b-X\c-X) J Y\j (b-Y)(c-Y) 

where \ x - \ Y = -{a + 2)EA+A-/16, a x = (a - 2)EX X /16, a Y = (a - 2)EX Y /1Q, 
b= -A+/A-, c= -A-/A+. 

A further change of coordinates 

X = - ysn 2 (a' + iK',k), Y = —ysn 2 (f3' + %K\ k) , k = ^± 

is convenient for writing the Schrodinger equation 



d 2 V 



Em 



(a + 2)k± (sn 2 (a' 5 k) - sn 2 ((3\ k)) + k 2 (a - 2) 

The separated equations are versions of Lame's equation ^2]- Indeed if we look for 
solutions of the form \l/ = A(a')B(/3') then 

+ (~^E{a + 2)snV, k) - A{ct) = , 

d 2 B{(3') + f_J_ k 4 E ^ a + 2) sn 2 (p', k) - \ 2 \ B(P') = 



d(3' z V 16 
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where Ai — A 2 = —E(a — 2)k 2 /16. Solutions of these separation equations can be repre- 
sented as Riemann P functions ^3] of the form 

/ k~ 2 oo \ 



P(z) 



1 (1 - \^A + k 2 E{a + 2)) sn 2 (z,k) 



i i 



1 l 

2 4 



l + \^A + k 2 E{a + 2) 



J 



\ 2 2 

for z = a', (3' . 

The case 7 = can easily be deduced by putting A + = iA_, as can be seen from the 
expression for the associated classical constant. 

If 7 = 2 then a convenient choice of coordinates is 



x = log^tan(<// — iu'yj , y = logftan(</?' + iu>') 
The corresponding classical Hamiltonian has the form 



H 



a + 2 a - 2 
+ 



The classical constant is 



X 1 +X 2 + 2K 2 = aH + 



sinh 2 2u>' sin 2 2ip' 



(a + 2) sin 2 2^'p 2 - (a - 2) sinh 2 2a/p 2 , 



(a + 2) sin 2 2ipf + (a - 2) sinh 2 2a/ 
The Hamilton- Jacobi equation in these coordinates is 



OS 
3y7 



fdS\ 



a+2 a -2 
+ 



E 



sinh 2a/ sin 2 2(p' 



which has solutions 



S(<p', J 



arctan 



(a - 2)E 
A 



sec 2 2ip' + tan 2 2ip' 



■ a/(o — 2)E arctantu / sec 2 2ip' + 



A 



[a - 2)E 



tan 2 2ip' 



ir- (a + 2)E , 2 l2 

H — v A arctan \ / sech 2uj' — tanh 2o/ 

2 V A 



; a/(o + 2)E arctanhi / sech 2 2u/ 



A 



{a + 2)E 



tanh 2 2oo' 



The corresponding Schrodinger equation is 

d 2 $> d 2 ^ 
+ 



dip 



duo 



i2 



a + 2 a 



Em 



sinh 2 2a/ sin 2 2(p' 



(4.7) 
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which has solutions of the form 



= Vsin2^sinh2^' Pp (a ~ ] (cos2^)^ V 1 j (cosh2^) 

where Py{z) is a solution of Legendre's equation. 

This completes the list of possible coordinate systems which are inequivalent and 
separable for this particular Hamiltonian. We notice in particular that the equation 
H — E = can be written in the equivalent forms 

Apl+Pl) + \E^a-2 + (a + 2)^j = 0, J 2 + J 2 + J 2 - E ^ J^^ + a^j = 0, 

both superintegrable systems on the complex two-sphere, the first of which is written in 
horospherical coordinates. 

4.2 Superintegrability for Darboux spaces of type four. 

There are various possibilities for the potential in this case: [A], [B], [C], [D]. 

[A1 " = ~ (a + 2)^a-2)* + °> + " 2 + h) + + * 

The additional constants of the motion have the form 
Rx = J fsT 2 + a 1 (/i 2 + z/ 2 ) + a 3 (/i 2 + z/ 2 ) 2 , 

2ai ((a + 2)/i 2 - (a - 2)z/ 2 ) + 16a 2 + 4a 3 ((a + 2)/i 4 - (a - 2)z/ 



i? 2 — -Aa + 



(a + 2)/i 2 + (a - 2)z/ 2 
The corresponding quadratic algebra relations are determined by 

R 2 = WR^l - 256a 3J R 2 - 64aii? 1 i? 2 - 256aa 3 HR 1 - 1024a 2 a 3 .Ri + 6Aa 1 HR 2 
- 256a 3 H 2 - 64ai(a + 2)H - 256a 2 a 2 . 

This Hamiltonian admits a separation of variables in coordinates corresponding to the 
equivalence first, second and third classes of Section l4~Tl For the second this is covered 
by the choice of coordinates //, v. 

(i) For coordinates corresponding to the first equivalence class, we obtain the Hamil- 
tonian in the form 

sin 2 2u (pi +p 2 v + 4aie 2l> + 4a 2 cosec 2 2w + Aa 3 e 4v ) 



H 



2 cos 2u + a 



(n) For coordinates corresponding to the third representative (|4.5jl the Hamiltonian 
takes form 

4(p 2 + p 2 ) + 4aic 2 (cosh 2 u — cos 2 ip) 



H 



+ 



- 2) (sech 2 u; — sec 2 cp) — (a + 2) (cosech 2 kj + cosec 2 <^) 
16a 2 (cosech 2 2a; + cosec 2 2y9) + a 3 c 4 (sinh 2 2u + sin 2 2(p) 
[a — 2) (sech 2 a; — sec 2 (p) — (a + 2) (cosech 2 ^ + cosec 2 ^) 
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The quantum algebra relations are 

[R,^] = -8{R U R 2 }- 16R 2 -32a x F, 

[R, R 2 ] = 8R 2 2 - 2hQa 3 R l - 128aa 3 H - 32(a 2 + 4a 3 + 16a 2 a 3 ) 

together with the operator relation 

R 2 = 8{R U R 2 2 } - 256a 3 i? 2 - 80R 2 2 - 256aa 3 HR 1 - 64(16a 2 a 3 + a\ + 4a 3 )i? 1 

+ 64 ai HR 2 - 25Qa 3 H 2 + 64a(4a 3 - a?)# + 128(al 2 + 4a 3 + 8a 2 a 3 - 2a 2 a 2 ). 

As in the case of free motion we observe that the equation H = E is 

2^2^ | a2-\(a-2)E a 2 -\(a + 2)E 
Pt + Pu + a i + — 2 + — 2 + W + v = » 

a superintegrable system in flat space with rearranged constants, that separates in elliptic 
and hyperbolic coordinates. 



[B] H 



sin 2u[p v +p u + . 2 + — -2- + 6i 
smh f cosh t> 



2 cos 2m + a 
The additional constants are 

Rl = Al + X 2 



26 x cosh 2w + (b 2 + b 3 ) (4 - a 2 ) + (cos 4m + 2a cos 2m + 3) 1 ? '~ 



+ 



R 2 = K 2 + — V + 



sinh 2 u cosh 2 t> 



2 cos 2u + a 



sinh 2 f cosh 2 v 
The quadratic algebra is given by 

R 2 = IQR\R 2 - UR\ - UaHR\ + 64(26 3 - 26 2 - b^Rj + 32a(b 2 + b 3 )R x R 2 

- 6AH 2 R 2 + 64(6 2 + b 3 )HR 1 + 128a(b 3 - b 2 )HR 2 

- 16 ((4 - a 2 )(b 2 + b 3 f + 8&i(6 2 - 6 3 ))i? 2 + 128(6 3 - h)H 2 - QAh(b 2 + b 3 f . 

This Hamiltonian admits a separation of variables in coordinate systems corresponding 
to the first and fourth equivalence classes of ()4.2|) . The defining expressions have already 
been given in terms of coordinates for the first. For the fourth, we distinguish two cases. 

(i) 7 ^ 2 

X{A.X - A + )(A + X - A_)p 2 x + Y(A„Y + A+){A+Y + AJ) V \ 



H = 16XY 



A + A_{X -Y)(a-2-(a+ 2)17) 



Ab 2 (A 2 -AV)XY Ab 3 (A 2 -A 2 ,)XY 

b 1 (XY + l)+ . . -p ~ ... +J V n l " 



+ ^_)(a + x + a„) (A_y + A+)(A_x + ^ 

a - 2 - (a + 2)XY 
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(n) 7 



, / x 1 \ 46 2 46 ? 



sinh 2u/ sin 2 2</// cos 2 2<// cosh 2a;' 
a + 2 a - 2 



+ 



sinh 2 2u/ sin 2 2^' 
The corresponding quantum algebra relations are 

[R, Ri] = -8Rj + 96R 2 2 + UaHR 2 - 16a(6 2 + b^R-y + 64(26 2 - 2b 3 + b x + 3)i? 2 
+ 32if 2 + 32a(26 2 - 2b 3 + l)H 

+ 646x(6 2 - 6 3 ) - 8(a 2 - 4)(6 2 + 6 3 ) 2 + 32(&x + 26 2 - 2o 3 ) , 
[R,R 2 ] = 8{R 1} R 2 } + 16a(b 2 + b 3 )R 2 -lQR 1 + 32(b 2 + b 3 )H -16a(b 2 + b 3 ), 

R 2 = -64R 3 2 + 8{R 2 , R 2 } - UaHR\ - 6AH 2 R 2 - 80R 2 - 64(26 2 - 2b 3 + 6 X + 7)i? 2 
+ 16a(6 2 + b 3 ){R u R 2 } + 64(6 2 + & 2 )M + 64a(2o 3 - 26 2 - 1)HR 2 

- 160a(6 2 + 63)^1 + 16 ((a 2 - 4)(6 2 + b 3 ) 2 + 8( 0l + l)(b 3 - b 2 ) - Ab x + 32^ R 

+ 128(63 - b 2 + l)H 2 + 128a(6 2 - b 3 + 1)# 

+ (&2 + &s) 2 (128 - 80a 2 - 646!) - 128(&i + 2)(6 3 - 6 2 - 1) - 256. 

As in the case of free motion, equation H — E = is 

-,9 ,n T o 26l 26 2 

Jl + J 2 2 + Jl + 7 r- + - r- + 6 3 

V4 + 4 (si + V4 + 4) V4 + 4 («i - V4 + 4) 

-e (** = + a )=o, 
\V4T4 J 

a superintegrable system on the complex sphere that separates variables in spherical, 
elliptic and degenerate elliptic type 1 coordinates. 

2 2 C l C 2 ( 1 1 

P 9 > + Pu> + — + 72— " + C 3 — — ~ . 2 

„ r cos z ip' cosh a;' \ sm V sinh a;' 
[CJ ii = 



2 



a + 2 a -2 
+ 



sinh 2u/ sin 2 2<p' 

These are coordinates associated with 7 = 2 in the fourth representative from (J4.2)) . The 
constants of the motion associated with this Hamiltonian are 

i?i = X 1 + X 2 + 2K 2 + aH 

a + 2 f c 3 Ci \ a — 2 ( c 3 c 2 



sinh 2 2uj' ysin 2 ^' cos 2 ip' J sin 2 2<p' \sinh 2 u/ cosh 2 u/ 

H a~+2 a^2 

+ 



sinh 2u>' sin u' 
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x l - x 2 + 



a + 2 a-2 



sinh 2 2uj' sin 2 uj' 



x 



+ 



-r—o I c i cos h 2o/ tan 2 ip' — c 2 cos 2ip' 



sinh 2uj' 
a-2 



c 3 ( 2 cos 2 <//(sinh 2 a/ - sin 2 cp') ) + 1 
sin 2 </?' 



sin 2^7 



- I c 2 cos 2yj' tanh uj' + Ci cosh 2a;' — 



C3 ( 2 cosh 2 ci/(sinh 2 uj' — sin 2 y?) + 1 



sinh uj' 



They satisfy the quadratic algebra determined by the identity 

R 2 = 16Rl - 16i?ii?l - 32aHRl + 32(c 2 - c^Rj + 16(a 2 - 4)if 2 i? 1 
+ 32((a + 2)ci - (a - 2)c 2 + Ac^HRt 

- 16(2c 2 - c 2 - c 2 + 6c 3 (c! + c 2 ) + 4cic 2 )i?i - 32(c 2 - c 3 )( Cl - 03)^2 

- 16((a + 2)(ci - c 3 ) 2 + (a - 2)(c 2 - c 3 ) 2 )# - 32( Cl - c 2 )(3c 2 - Cl c 2 - c 3 (ci + c 2 )) 

The Hamiltonian admits a separation of variables in a number of coordinates systems 
corresponding to various combinations of the operators R\ and i? 2 . We exhibit the various 
possibilities. 

(i) For the constant R\ — i? 2 , the associated separable coordinates are those corre- 
sponding to the third representative in ()4.2|) with 7=1. In these coordinates, the 
Hamiltonian is 



H 



Ci + c 2 + 2c 3 (ci + c 2 ) cosh 2uj c 3 cos 2ip 



2 sinh 2 2uj 



+ 



sin 2 93 



(a-2) 



1 



cosh c<j cos 2 ip 



(a + 2) 



sinh sin ip 



(ii) In coordinates corresponding to rotations of the fourth representative in ()4.2)1 with 
7 ^ 0, 2, that is, B\X X + {B\ - B 2 )X 2 + {2B\ - B 2 _)K 2 , the corresponding Hamil- 
tonian has the form 



H = 16 



+ X)(B ± + X)p| + Y(B T + y)(B± + 

1 1 1 



ci / 1 1 



C 3, 



2 B 4 2 



+ 



4 \r x 
1 



(Bl - B\ 



1 + B T Y 1 + B T X 1 + B±Y 1 + B±X 
a-2 a-2 a + 2 a + 2 



+ 



l + B±X l + B ± Y 1 + B T Y l + B^X 
a — 2 a — 2 

+ (a + 2)(X-Y) 



X 



Y 



Here, B± = B + /B_ and B T = B-/B+. The Hamiltonian associated with i? 2 can 
be obtained from this last case by taking j3_ = \/2B + . 
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The quantum algebra relations are 

[R, Ri] = 8{Ei,E 2 } + 16E 2 + 16(ci-c 3 )(c 2 -c 3 ), 
[R,R 2 ] = 24R 2 1 -8R 2 2 - 32aHR 1 + 8(a 2 -A)H 2 + 32( Cl - c 2 - 
+ 16((a + 2)ci - (a - 2)c 2 + a + 64c 3 )E 
+ 8c 2 + 8c 2 - 16c 2 - 32cic 2 - 48c 3 ( Cl + c 2 ) + 16(ci - c 2 ) . 

The operator identity is 

E 2 = 16E? - 8{R U E 2 } + 32(c 2 -c x - |)E 2 - 80E 2 . + 16(a 2 - 4)# 2 Ei 

+ 32 ((a + 2)ci - (a - 2)c 2 + 4c 3 + a) Mi + 16 (c 2 + c 2 - 2c 2 ; - 6c 3 (ci + c 2 ) 

- 4 Cl c 2 + 2(ci - c 2 ) - 8) - 32(c 2 - c 3 )(ci - c 3 )E 2 + 16(a 2 - 4)E 2 

- 16 ((a + 2)((ci - c 3 ) 2 - 2ci) + (a - 2)((c 2 - c 3 ) 2 + 2c 2 ) - 8c 3 - AajH 

- 32(ci - c 2 )(3c 2 - cic 2 - c 3 (ci + c 2 )) 

+ 32(c 2 + c 2 , - 4c 3 (ci + c 2 ) - 2cic 2 + 2ci - 2c 2 ). 

As in the case of free motion, the equation H = E is 

r2 , t2 , r2 i(ci + c 2 + 2c 3 )si i(ci - c 2 )(si + is 2 - s 3 ) 

Jl + J 2 + J 3 ^ ^ + 



4 V / S 2 r Tsf 4^^(51 + is 2)(Ss -«S2y 

(2c 3 - ci - c 2 )(si + is 2 + gg) , i(ci - c 2 ) / 2igi 

H = 1 -= hi a H ^=^= |=U. 

y / (s 1 + is 2 )(s 3 + is 2 ) 4V2 y V s 2 + S i, 

which is a superintegrable system on the complex sphere, with rearranged constants, that 
separates variables in elliptic and degenerate elliptic coordinates of type 1. 



4//V 



[D] H = 



(a + 2)/i 2 + (a - 2)u 2 
This Hamiltonian admits three classical constants of the motion 

d(fi 2 + is 2 ) 2 „ v 16d 

Ei = Xi + - — -— , R 2 = X 2 + - — — -J , K = fj l pt l + vp v . 

(a + 2)[i 2 + (a — 2)v 2 (a + 2)fj J 2 + (a — 2)u 2 

The Poisson quadratic algebra satisfies the relations 

{K,R 1 } = 2R 1 , {E,E 2 } = -2E 2 , {R u R 2 } = -8K 3 - AaKH - lQdK. 

These three extra constants are related via the identity 

—R 1 R 2 + E 4 + aHK 2 + UK 2 + H 2 = . 

This Hamiltonian admits a separation of variables in all the coordinate systems that 
are possible. We need only give the expressions in terms of the fourth representatives. In 
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the coordinate system associated with the fourth representative and for which 7 7^ 2 the 
Hamiltonian can be written as 

x(A + x - a_)(a_x - A + )p\ - + + a + ) p 2 y 

(X -Y)(a-2-XY(a + 2))A + A^ 
AdA + A_ (X 2 Y + Y + XY 2 + X) 



(X -Y)(-a + 2 + XY(a + 2))A + A^ 
and for the case 7 = 2 this Hamiltonian has the form 



Pt' +Pt> + d[ . -2- - + 



H 



sinh 2u' sin 2 2(p' 



a+2 a-2 



+ 



sinh 2 2tu' sin 2 2<p' 
The corresponding quadratic algebra relations are 

[K,R 1 ]=2R 1 , [K,R 2 ] = -2R 2 , [R h R 2 ] = -8K 3 - AaHK - IQdK - 4K, 
subject to the operator identity 

--{Ri, R 2 } + H 2 + aHK 2 + K 4 + aH + (5 + Ad)K 2 + Ad = . 

This completes the analysis of the superintegrable potentials associated with the four 
metrics of Darboux. 



5 Relationship to constant curvature superintegrable 
potentials 

In sections EH3] we have found, by means of exhaustive calculation, all superintegrable 
potentials in the Darboux spaces of revolution having two or more quadratic integrals. 
Once these are expressed in suitable coordinates, it is clear that each is simply a multiple 
of one of the superintegrable potentials on the complex Euclidean plane or two-sphere, 
that have been enumerated in [7], though that was by no means evident in advance. 
In each case we can start with a Hamiltonian of the form 

H = H + aV , (5.1) 

where Vq is a function of the coordinates x and y, and a is a constant. Dividing 
the Hamilton- Jacobi equation, H = E, throughout by Vq and rearranging gives a new 
Hamilton- Jacobi equation in which the roles of the energy E and parameter a have been 
exchanged. 

Clearly, the integrability and separability of one system guarantees that of the other. It 
is this relationship between the harmonic oscillator potential written in Cartesian coor- 
dinates and the Coulomb potential in parabolic coordinates that has been discovered by 
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many authors. Transformations of this type relating integrable systems were described in 
a more general context by Hietarinta et al in and called coupling constant metamor- 
phosis. See also ^2] where the Stdckel transform and its close connection with variable 
separation was emphasized. 

The preservation of integrability under such a transformation can be demonstrated 
explicitly by noting that if {H , L } = and 

H = H + aV and L = L + a£ (5-3) 

are in involution, i.e., {H, L} = 0, then so are 

H'=^y- and L' = L — i H' . (5.4) 

Any identities involving integrals associated with ([5.1)1 . give rise to corresponding 
identities involving integrals associated with (|5.2j) and are obtained by the replacements 

a -> -H' and H -> . (5.5) 

5.1 Generating the Darboux spaces of revolution by coupling 
constant metamorphosis 

Taking each of the degenerate potentials from [Jj, that is, the potentials with Hamiltonians 
having one first order and two quadratic integrals and performing a coupling constant 
metamorphosis we arrive at a Hamiltonian having one first order K and two quadratic 
constants, X\ and X 2 . These must be free Hamiltonians either on one of the four Darboux 
spaces of revolution or one of the constant curvature spaces, E 2 (C) or ^(C). After 
comparing the Hamiltonians so generated, it can been seen that this approach generates 
all of the Darboux spaces of revolution. 

Knowing the Poisson algebra for each Hamiltonian involved and how coupling con- 
stant metamorphosis modifies this algebra, we can determine which Hamiltonian has been 
generated, even if it appears in unfamiliar coordinates. Note that some transformations 
reproduce the free Hamiltonian on E 2 (C) or S 2 ,c, and some Darboux spaces can be gen- 
erated from two distinct constant curvature potentials. 

For each Hamiltonian we have four linearly independent constants of the motion. 
These, however, cannot be functionally independent and there is always a polynomial 
identity in K, X±, X 2 and H that is of fourth order in the momenta. We can use this 
identity to classify the possible Hamiltonians. Up to freedoms in choosing X% and X 2 , 
scalings of K and coupling constant metamorphosis, we find that there are 5 classes of 
identities that involve all of the constants. The correspondences between these identities, 
degenerate superintegrable potentials from [7j and the Darboux spaces of revolution are 
summarised in Table Note that because we allow coupling constant metamorphosis, 
H has the same status as parameters in the the potential and the coefficients A and B 
appearing in the representative identities may be functions of H. The labels in bold (e.g. 
E3, S3,. . . ) refer to 7J. Those Hamiltonians in Tabled on the complex two-sphere, that 
is, S3, S5 and S6, are represented with three coordinates Si, s 2 and s 3 constrained by 
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4 + 4 + 4 = 1 and Ji = s 2Ps 3 - s 3 p S2 , J2 = s 3 p Sl - s 1 p s:j and J 3 = s x p S2 - s 2 p si . The 

4+P 2 y 



potentials E12, E14, E4 and E13 are functions of x — iy and hence division of p 2 + pi 



by these potentials reproduces the flat space Hamiltonian. 

For example, starting from the algebraic identity for constants associated with the 
Hamiltonian and integrals 

H = p 2 x +p 2 y + ^ + a (E6), (5.6) 

Xi = [xpy - yp x )p x - ^| , X 2 = (xp y - yp x ) 2 + ^~, K = p y , 
that is [Zj, 

X\ + K 2 X 2 — (H — a)X 2 + aK 2 = , 
we find that applying the transformation (|5.4jl gives 

H = J_ + i ' X i = ( x Py ~ VPx)Px ~ -5 j_ ~ x = +xp x p y , 



X' 2 = (xpy - yp x ) 2 + — \ y = 5—^ - 2xyp x p y , K' = K , 

-'- -g- + 1 X 1 + 1 

Xf + X /2 X 2 - H'X 2 - H'K' 2 = . 



x 2 \ + 1 + 1 

as* 

and using (|5.5j) . 



Then 

X'[ = 2X[ , X' 2 ' = -X' 2 + , H" = H' , K" = K' 

gives 

X" 2 - 4X" 2 X2 + AH"X 2 - AH" 2 = , 
the identity (|2.3|) associated with the Darboux space of type two (|2.1j) . 



5.2 Generating superintegrable potentials on Darboux spaces 

The in equation ()5.3)1 may itself contain potential terms and if these are chosen so 
that H is superintegrable, then so will be H' . 

For example, taking the superintegrable Hamiltonian on the complex two-sphere SI 

tt 72, t2 , ft, 01 , l 3 ^ 7(1 

-H = J x + J 2 + ^3 + 7 : — 77 + 7 : — 77 + 7 : — 77 + ° 

(S 1 -IS 2 ) 2 (S!-ZS 2 ) 3 (Sx-^Ss) 4 

and dividing though by (s 1 — is 2 )~ 2 — 1 gives, after a change of coordinates, the super- 
integrable potential [A] in a Darboux space of type 2. The same Hamiltonian can be 
generated by dividing E2 throughout by x~ 2 + 1. 

Each potential in Table d is compatible with the addition of further terms while 
maintaining superintegrability, and in using the method demonstrated above, all superin- 
tegrable Hamiltonians found in Sections I2HU can be generated. The correspondences are 
given below. 
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Table 1: Correspondences between constant curvature superintegrable potentials and 
Hamiltonians for Darboux spaces of revolution. 



Degenerate superintegrable 
potential on -E^C) or 
S 2 (C) 


Hamiltonian for Darboux 
space of revolution 


Representative identity 


E5: 4x 


n . PI+PI 
4u 


X 2 + AX 2 + X 4 + B = 


E6: -1 + 1 

x A 

S5: 

(si - is 2 y 


u\pl + p 2 v ) 

u 2 + 1 


X 2 + K 2 X 2 + AX 2 + B = 


E12: — -i ^ + 13 

y{x — iy) 2 + c 2 

E14: - a , + /3 


E 2 (C) 


X 2 + K 2 X 2 + A = 


E3: + + 4 
E18: — — + 1 

V x2 + y 2 


3 ' 4 + W 2 +W 2 


X X X 2 + AX 2 + 5 = 


S3: a4 ^ 2 a + 2 

4 

bo: — ; h a 

V S l + S 2 


p. 

4 ' a+2 , a-2 
9 — r~ 9 


XxX 2 + K 4 + AK 2 + 5 = 


E4: a(x - iy) + (3 
E13: , ^ , + /3 


£ 2 (C) 


K 2 X 1 + AX 2 + 5 = 



38 



5.2.1 Darboux spaces of type 1 

The potential E5, V$ = Ax, appears in each of 

E2 : a(Ax 2 + y 2 )+(3x+^ + 5 

x 2 

E3' : a(x 2 + y 2 ) + (3x + yy + 5 

Oi 7(2x — iy) 

E9 : , +(3x+ ! \ y ' - + 5. 

Vx - iy yjx — iy 

The potential labelled E3' is a translation of E3 . Adding these potentials to H = pl+p 2 
and dividing by 4x produces the two real non-degenerate potentials found in pQ and an 
additional complex one given in this paper. (The details of the quadratic algebra and 
defining operators for the Hamiltonian derived from E9 can be computed using (|5.2jl ). 

5.2.2 Darboux spaces of type 2 

The potentials E6 and S5 appear in each of following. 

El [B] : a (x 2 + y 2 ) + ^ + ^ + 6 



x- r 



7 



E2 [A] : a(Ax 2 +y 2 )+(3x + ^ + 5 

yA 

E16 [C] : . 1 (a+ ^ + } ) ■■■<) 

a/x 2 + y 2 \ x + yx 2 + y 2 x — a/x 2 + y 2 

51 [A] • a | 1333 | ^( 1 ~ 4 ^) | S 

(si-is 2 ) 2 (si-is 2 ) 3 {x-iyY 

52 [B] : « + P^ + p±^ + S 

S4 [°] : ( a - ^ + -7^=f + 7 . 7 r2 — 2 + 5 

The superintegrable system generated after dividing by x~ 2 + 1 or (si — is 2 )~ 2 — 1 as 
appropriate is indicated by label the [A], [B] or [C]. The apparent over abundance of 
superintegrable potentials generated in this way for D 2 is resolved by noting that the 
same potential can appear in more than one coordinate system. 

5.2.3 Darboux spaces of type 3 

The potentials E3 and E18 appear in each of 

El [B] : a (x 2 + y 2 ) + ^ + ^ + 5 

x y 

E3' [A] : a{x 2 + y 2 ) + f3x + ^y + 5 

a(x — iy) (3(x + iy) 



E7 [D] 



^ a/(x - iy) 2 - c 2 ((x - iy) + \/{x - iy) 
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+7(a; 2 + y 2 ) + 5 
-r^ r^, n a(x + iy) 8 , 9 9x 

E16 [B] : 1 |g+ ^ + I ) 

y x 2 + t/ 2 \ a; + a/x 2 + y 2 x — a/x 2 + y 2 

E17 [C] : . " + , P . + T , + 5 

a/x 2 + y 2 (x + n/) 2 (x + ^)A/x 2 + y 2 

E19 [D] : v u ' = + 



a/(x-u/) 2 -4 yj (x + iy){x - iy + 2) 

+— f 7 ^= + £ 

V / (ar + «j/)(z-ij/-2) 

E20 [A] : 1 (g + /3A/x+ yG 2 + y 2 + -yJ x - ^/x 2 + y 2 ) +(5 
a/x 2 + y 2 V / 

As before, once the possibility of changes of coordinates in taken into account, the above 
list produces only those superintegrable potentials found in section 13.21 

5.2.4 Darboux spaces of type 4 

The potentials S3 and S6 appear in each of 

S4 [A] : 7-^ + ^^ + ^-^= + 6 

{si-ts 2 ) 2 yjs\ + s 2 2 (s l -is 2 )yjs\ + s 2 2 

57 [B,C] : -^^ + ^ +^ 2+ 5 

58 [C] • aSl + + is 2 - gg) + 7(fi + ^2 + s 3 ) + tf 

a/s| + s§ a/(si + is 2 )(s 3 -is 2 ) \/{si + is 2 ){s 3 + is 2 ) 

59 [B] : ^ + | + J + 5 

6 1 5 2 5 3 

As before, once the possibility of changes of coordinates in taken into account, the above 
list produces only those superintegrable potentials found in section 14.21 



6 Conclusion 

In this paper we have discussed in some detail three of the four Darboux spaces of rev- 
olution that have at least two integrals of classical motion quadratic in the momenta in 
addition to the Hamiltonian. In each case we have also presented an exhaustive list of 
potentials for each of these spaces which when added to the Hamiltonians of these spaces 
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preserve this property i.e. that there are still two extra integrals of the classical motion. 
These are the superintegrable systems associated with the systems of Darboux. The 
property of extra integrals also extends easily to the case of the corresponding quantum 
systems. For each of these systems we have calculated the corresponding quadratic algebra 
relations and shown that in each case the Hamiltonians that we obtain arise from constant 
curvature systems via a coupling constant transformation. We have also discussed the 
solutions of the corresponding classical and quantum problems in each of the inequivalent 
coordinate systems and have also given some of the embeddings of these spaces in three 
dimensions. In the last section we have shown how the free Hamiltonians of Darboux are 
related to the superintegrable Hamiltonians on spaces of constant curvature via coupling 
constant transformations. We also list how the corresponding superintegrable systems of 
spaces of constant curvature are related in this way to the superintegrable systems that 
we have found. This classification is comprehensive and complete. 

Let us very briefly review the current status of superintegrability in two-dimensional 
spaces. Most of the published work [SI IH El El El concerns quadratic superintegrability 
for classical, or quantum Hamiltonians of the form kinetic energy plus a scalar potential. 
Once a specific space is chosen, superintegrable systems in the space can be classified 
under the action of the corresponding isometry group. Systems in the same class are not 
only mathematically equivalent, but also have the same physical properties. In classical 
mechanics they will have the same trajectories and the trajectories will be periodic, if 
they are bounded. Similarly, in quantum mechanics superintegrable systems in the same 
class will have the same energy levels and eigenspaces. 

Quadratically superintegrable systems exist in spaces of constant curvature and also 
in Darboux spaces. A Darboux space is defined by the fact that it allows one Killing 
vector and two (irreducible) Killing tensors. This article completes the task of classifying 
all quadratically superintegrable systems in all of the above spaces. 

The results are quite rich. Indeed, in the real Euclidean space E 2 , we have four 
£7(2) classes of superintegrable systems [HI El- They are physically quite diverse. One is 
an isotropic harmonic oscillator with additional terms, called El above in Section El A 
second is an anisotropic harmonic oscillator with additional terms (called E2 above). The 
third and fourth are Kepler (or Coulomb) systems with two different types of additional 
terms, respectively. In complex Euclidean space £7 2 (C), or correspondingly in the pseudo- 
Euclidean space £7(1, 1), one obtains 6 more classes [3]. 

Two classes of superintegrable systems exist on the real sphere S 2 , four more on the 
complex sphere S 2 (C) [Hj. On the real Darboux spaces D\, . . . , D4 we have obtained 3, 4, 
4, and 4 classes of systems, respectively. One more for the complex space .03(C). 

From the mathematical point of view the situation is much more unified. As was 
stressed above, superintegrable systems that may correspond to quite different physical 
situations may be related by coupling constant metamorphosis. Once we allow this type of 
equivalence, many fewer equivalence classes exist. For instance, in real Euclidean space we 
only have two classes, because the Kepler potentials with additional terms are equivalent 
to isotropic harmonic oscillators (in one case with the additional terms). All superin- 
tegrable systems in Darboux spaces are related by coupling constant metamorphosis to 
systems in spaces of constant curvature. For D\, D 2 and D 3 this is always flat space, 
complex or real. Two of the systems in D4 are related to systems in real Euclidean space. 
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The other two are related to systems on a complex sphere. The relation is of course not 
unique and depends on the choice of coordinates (see Sectional). 

A typical feature of quadratic superintegrability for scalar potentials is that quantum 
and classical superintegrable potentials coincide. They allow separation of variables in at 
least two coordinate systems in the Schrodinger and Hamilton- Jacobi equation, respec- 
tively. 

Superintegrability involving third order integrals of motion has also been considered 
|16| ITT]. There the situation is quite different. Multiseparability is lost. More interestingly, 
quantum superintegrable systems exist (in real Euclidean space) that have no classical 
analog (in the classical limit they reduce to free motion). 
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